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ABSTRACT 

We present the detection and analysis of line index variations in the pulsating sdB star 
PG 1605+072. We have found a strong dependence of line index amplitude on Balmer 
line order, with high-order Balmer line amplitudes up to 10 times larger than H/3. 
Using a simple model, we have found that the line index may not only be dependent 
on temperature, as is usually assumed for oscillating stars, but also on surface gravity. 
This information will provide another set of observables that can be used for mode 
identification of sdBs. 
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1 INTRODUCTION 



Subdwarf B (sdB) stars are hot, core-helium burning 
stars, with hydrogen envelopes that are too thin to sus- 
tain nuclear fusion. The discovery of multimode pulsa- 
tions in some sdBs should allow the use of asteroseis- 
mology to probe their atmospheres, and thereby help 
to answer the many questions remaining about sdB for- 
mation and evolution. 

Until recently, time-series observations of pulsating 
sdBs have been limited to photometry. We have shown 
that it is possible to detect radial velocity variations 
of Balmer lines in the pulsating sdB star PG 1605+072 
using 2m-class telescopes (O'Toole et al. 2000, 2002, 
hereafter Papers I and II). Four-metre class telescopes 
have also been used to detect velocity variations in sdBs 
(Jeffery & Pollacco 2000; Woolf et al. 2002). In Paper 
II we found closely spaced peaks in the amplitude spec- 
trum of PG 1605+072, as well as amplitude variation 



* Based on observations made with the Danish 1.54 m telescope 
at ESO, La Silla, Chile. Part of the data presented here have 
been taken using ALFOSC, which is owned by the Instituto de 
Astrofisica de Andalucia (IAA) and operated at the Nordic Op- 
tical Telescope under agreement between IAA and the NBIfAFG 
of the Astronomical Observatory of Copenhagen. 



in at least one of these peaks. In this paper we de- 
scribe the analysis of line-index variations of Balmer 
lines in PG 1605+072. The atmospheric parameters of 
this sdB have been determined from high-resolution 
optical spectra and NLTE model atmospheres to be 
T G ff =32 300+300 K and log g=5. 25+0. 05 (Heber et al. 
1999). The helium abundance was found to be subsolar, 
i.e., log (He/H)=— 2.53+0.1. The star was also found 
to be rotating with a projected rotational velocity of 
Sgkms" 1 . 

The use of equivalent widths of Balmer lines to 
monitor stellar oscillations of solar-type stars was first 
proposed by Kjeldsen et al. (1995). The Balmer lines 
are very temperature sensitive - the equivalent width of 
these lines depends on the number of hydrogen atoms 
in the second level. According to Kurucz (1979), the 
maximum line strength of Ha to H5 is between 8500 K 
and 9000 K at log 5=4.0. Since the temperature changes 
slightly with the oscillations, the equivalent width must 
also change. The idea was extended by Bedding et al. 
(1996) to include metal lines. In this paper we calcu- 
late line indices, which depend on the temperature and 
surface gravity sensitivity of the equivalent widths of 
spectral lines. This method has been used to identify 
modes in the 8 Scuti stars FG Vir (Viskum et al. 1998) 
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Figure 1. Spectrum of PG 1605+072 taken using DFOSC at the 
Danish 1.54 m telescope. The continuum is clear until around He. 



and BN Cnc (Dall et al. 2002), and to detect equivalent- 
width oscillations in the Ha line of the roAp star a Cir 
(Baldry et al. 1999). Here we present the first detection 
of line index variations in a pulsating sdB star, namely 
PG 1605+072, and find a dramatic dependence of line- 
index amplitude on Balmcr line number. We show that 
a simple model can qualitatively explain this effect as a 
combination of temperature and surface gravity fluctu- 
ations. 



2 OBSERVATIONS AND REDUCTIONS 

In this paper we use the spectroscopic observations de- 
scribed in Papers I and II. Briefly, they consist of 16 
nights in May 2000 using DFOSC on the Danish 1.54m 
in Chile, 4 nights overlapping this using ALFOSC on the 
NOT 2.56 m and, to improve frequency resolution, a fur- 
ther ~1 hour per night for a month in March- April 2000 
using DFOSC. We also have 10 nights of observations 
using DFOSC from July 1999 using the same setup, as 
well as 2 nights using the coude spectrograph mounted 
on the Mount Stromlo 74-inch telescope in Australia. 
Finally, we have 12 nights of Johnson B photometry 
from the South African Astronomical Observatory, also 
taken in July 1999 (see Paper I). 

As described in Papers I and II, reductions of the 
spectra were done using standard IRAF routines for 
bias subtraction, flat fielding and background light sub- 
traction, and spectra were extracted using a variance 
weighting algorithm. A typical spectrum is shown in 
Figure 1. For details of photometric reductions see, for 
example, Koen et al. (1997). 
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Figure 2. Line index as a function of time for the H8 line. 

3 THE LINE INDEX 

We measure the line index, a quantity which depends 
on the temperature and surface gravity sensitivity of 
the equivalent width. It is defined as 

a l = E (i) 

(Dall et al. 2002), where C x is the continuum count level 
at pixel x, S x is the source count level, and W x is a 
gaussian-like weighting functionfor the line in question. 
The value of C x is found using C x = £ S X W% / £ W x 
where W x is a broad filter function and is also gaussian- 
like. W x is centered by moving an integration weight fil- 
ter across the line in question. The position that max- 
imises the sum through this filter is taken as the line 
centre (Dall et al. 2002). The width of the function can 
be chosen arbitrarily; we have chosen the approximate 
full width at half maximum of each Balmer line, which 
means we use a different filter for each line. 

Prior to calculating the line index, each spectrum 
was normalised to an average number of counts per pixel 
and divided by the local continuum. The line index was 
calculated by multiplying each Balmer line and the ad- 
jacent continuum by the weighting function (a supcr- 
super-gaussian: e~ x ), dividing the line by the average 
of the two continua, and then summing over the line. 
In this way, the line index can be thought of as similar 
to the Stromgren H/3 index, calculated using software 
filters. (For a justification of the term "line index" as 
an analogue to a colour index, see Dall et al. 2002). The 
software for calculating line indices has been collected 
in a software package called Ix, as described by Dall 
(2000). 

3.1 Individual Balmer Lines 

We have calculated the line indices of all Balmer lines 
from H/3 to H9. These values were normalised to a mean 
of zero and divided by the total mean to give the frac- 
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tional change (i.e. we calculated (A — (A))/(A), where 
(A) is the mean line index). A sample line index curve 
(of H8) is shown in Figure 2; the peak-to-peak scatter of 
the best quality data (from the NOT) is around 20%. 
As in Papers I and II, the quality of the data varies 
from night to night, so we have weighted the Fourier 
transform using the internal scatter. 

The amplitude spectra for the six Balmer lines from 
H/3 to H9 is shown in Figure 3. The left column shows re- 
sults for the 1999 data (Paper I), while the right column 
shows the 2000 data (Paper II). An amplitude change 
is evident over the 300 days between observations, sim- 
ilar to that seen in velocity. Also, line-index amplitudes 
vary strongly with Balmer line number in both sets of 
observations. We now consider possible reasons for the 
latter effect. 

3.2 Dependence of amplitude of Balmer-line 
order 

Could the variation with Balmer line order seen in Fig- 
ure 3 be an artifact of our measurement technique? 
There are several factors we should consider. 

Firstly, in sdB stars the Balmer lines are broadened 
by temperature and gravity, so it is difficult to measure 
the continuum near high-order lines (n=9-12), since it 
does not really exist in these regions. However, the am- 
plitudes of the strongest peaks differ by about a factor 
of three between H/3 and He, and these lines have good 
continuum, as seen in Figure 1. Still, the wings of the 
Balmer lines are very broad in sdBs and are dominated 
by gravity effects, so we will investigate this "pseudo- 
continuum" at high order Balmer lines in Section 4. 

A second explanation may be that the wings of the 
Balmer lines are oscillating out of phase with the cores. 
Since the wings extend out so far, they may be inadver- 
tantly sampled as part of the continuum and thereby 
boost the observed amplitude of line index change. We 
have determined the line index variation across each 
Balmer line at each of the frequencies in Table 7 of Pa- 
per II. The phases at each frequency are constant when 
we examine different slices across the line, leading us to 
rule out this explanation. 



4 A SIMPLE MODEL 

To investigate the observed behaviour, we have used H- 
Hc line blanketed, metal-free NLTE model sdB spec- 
tra (Napiwotzki 1997) at log.g=5.0, 5.25, 5.5 and 
T cff =30 000 K, 32 500 K and 35 000 K. The helium abun- 
dance is fixed at log(Hc/H)= —2.5. The calculations 
that follow were also carried out using metal line- 
blanketed LTE model sdB spectra with solar metallicity 
and Kurucz' ATLAS6 Opacity Distribution Functions 
(Heber et al. 2000), but since the basic conclusions are 
the same, we present only the NLTE results. A detailed 
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Figure 4. Model sdB spectrum with T of f=32 500K and 
log g=5.25, used to determine line index variations. The contin- 
uum is based on a typical NOT spectrum. Flux is in arbitrary 
units. 

discussion of NLTE vs. LTE models can be found in 
Heber et al. (2000). In the following sections we will 
compare the line indices of Balmer lines in these spec- 
tra with those of our observations. 

4.1 The Effect of Temperature 

For most oscillating stars, it is usually assumed that 
variations are predominantly effected by temperature 
changes, and that radius changes (and hence surface 
gravity changes) are negligible (e.g., Kjcldscn et al. 
1995). In this case, the oscillation amplitude in line in- 
dex is related to the temperature fluctuations by the 
relation 

SA _ OlogA ST 

A dlogT T ' U 

For cool stars 9 log A/9 log T is positive, but for hot 
stars such as sdBs, it is expected to be negative. (Note 
that for A stars it is approximately zero). To investigate 
the behaviour of the line index with varying tempera- 
ture for each Balmer line, we show in Figure 5 both A 
and 9 log A/9 log T as a functions of temperature. Here 
we have used LTE models with a larger range of ef- 
fective temperatures and the surface gravity fixed at 
log 5=5.5. Prior to calculating the line indices, we re- 
binned the model spectra to the approximate dispersion 
of our observed spectra and multiplied these spectra by 
a typical continuum from the NOT observations. We 
also used exactly the same parameters in our line index 
software for both models and observations. 

Figure 5 is similar to Figure 2 of Bedding et al. 
(1996), which was made for cooler main sequence stars. 
An important difference is the scale of change of line 
index with temperature. For cool main sequence stars 
there is a factor of three difference over around 1000 K 
for H/3, while for our sdB models, the difference for H/3 
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Figure 5. LTE model line index, A {top panel) and its deriva- 
tive d log A/<9 log T {bottom panel), both as a function of effective 
temperature. The symbols have the following meaning: diamonds 
- H/3; triangles — H7; squares — H<5; crosses - He; open circles - 
H8; stars - H9; filled circles - H10. 



Figure 6. line index of the highest peak (2742.85 /iHz) in the 
2000 amplitude spectrum compared with a simple model with 
temperature variations of -~500K (dashed line). 



is about a factor of two over 20 000 K. The line index of 
all Balmer lines except H9 and H10 seems to converge 
at about 40 000 K, although without modelling beyond 
this temperature we cannot be certain (and ~40 000K 
is the upper limit of the sdB regime anyway). As is 
expected, the change in line index with temperature is 
negative, that is, the line index of hotter sdBs is smaller. 

To investigate the "pseudo-continuum" idea sug- 
gested in Section 3.2, we rebinned our NLTE model 
spectra to the approximate dispersion of our observed 
spectra and again multiplied the resultant spectra by 
a typical continuum from the NOT observations. Our 
goal was to use a model that was as close to the ob- 
servations as possible. We calculated line indices for 
each of our models, once again using the same param- 
eters as for the observations. To simulate oscillations, 
wc took the difference between the line index of each 
model Balmer line at each of the different tempera- 
tures listed above. The surface gravity of the models 
was fixed at logg=5.5. The change in line index, AA, 
is approximately linearly proportional to the change 
in temperature AT c s . Because of this relationship, wc 
can scale the line index "amplitude" to a more realis- 
tic value of AT. We chose AT ~ 500 K, which is an 
order-of-magnitude estimate based on photometric am- 
plitudes. Finally, we normalised these differences by the 
line index values at 32 500 K (the approximate effective 
temperature of PG 1605+072), allowing us to compare 
fractional changes directly with our observations. 

The model line-index amplitude is plotted as the 
dashed line in Figure 6. Also shown is the line-index 
amplitude of the 2742.85 /iHz peak from the March- 
May 2000 observations. Both model and observations 
show an upward trend moving to bluer wavelengths, but 
the amplitudes do not match. A similar result is seen 
for each of the other frequencies in Table 7 of Paper II. 
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Figure 7. LTE model line index, A {top panel) and its deriva- 
tive d log A/<9 log g {bottom panel), both as a function of surface 
gravity. Symbols have the same meaning as Figure 5. 



Clearly then, these unusual variations depend on more 
than just temperature, particularly for the higher-order 
Balmer lines. We now consider the effect of changing 
surface gravity as well as temperature. 



4.2 Add Gravity and Stir 

Here we assume that the relationship between change 
in surface gravity and change in line index is similar to 
Equation 2, and that we can use the simple formula 



SA 

T 



6A\ 

"a; 



T c ff 



SA\ 



(3) 



logs 



to combine the effects of T e s and log g. The first term on 
the right of Equation 3 is the line index at fixed T c g, and 
the second term is at fixed logy. This linearity assumes 
that the oscillations are adiabatic, which implies that 
T e g and logy are in phase. 
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We show A and d log A/d log g as functions of log g 
in Figure 7. The important thing to note about this 
figure is that the line index of low-order Balmcr lines 
stays roughly constant with log<?, while for the higher 
order lines the line index decreases dramatically. This is 
related to the psuedo-continuum effect discussed earlier. 
The implication of Figures 5 and 7 is that line index 
amplitudes should be highest for higher order Balmer 
lines, which is what we see. We can now try to quantify 
this effect. 

Using a similar method to that used for tempera- 
ture changes above, but keeping T e ff fixed at 32 500 K, 
we have calculated the approximate effect of changing 
log g. To fit the combination of these effects, we min- 
imised the rms scatter of the difference between obser- 
vations and the model. The best-fitting combination of 
the temperature and surface gravity effects is shown in 
Figure 8. The fit appears to be, in general, very good, 
although there is some discrepancy at H10. For this fre- 
quency, we find that AT - 372 K and Alog.g - 0.047. 
Whether or not these values are reasonable will be dis- 
cussed in the next section. 

Assuming that this method is valid, we have deter- 
mined the temperature and surface gravity changes for 
each frequency in Table 6 of Paper II, as well as the 
frequencies measured from the 1999 observations (Ta- 
ble 5 of Paper II), and these are shown in Table 1. The 
velocity amplitudes derived in Paper II are shown in 
the last column. The uncertainties in both AT c ff and 
Alogg were found by fixing one of the (<5A/A)j, and 
changing the other until the rms increased by the mean 
error of the measured line index amplitudes. We have 
redone our calculations including the effects of rota- 
tion (vsmi = Sgkms" 1 for PG 1605+072; Heber et al. 
1999), and find that the values in Table 1 remain the 
same within 0.1%. 

In Figure 9 we show AT as a function of Alogg 
for both sets of observations, with la error ellipses de- 
termined from linear regression. Most of the points lie 
on or near the straight line. The slope of the line is 
8471 Kdex -1 . This can be compared with future mod- 
els. If there is linearity, it would be consistent with all 
modes having the same I value. Modelling by Kawaler 
(1999) found 4 modes to be I — 1 and one to be 
I = 2. We have detected 4 of the 5 modes predicted 
by Kawaler, including the mode predicted to be I = 2 
(around 2102 /iHz). In Figure 9, this mode appears to 
lie on the same line as the other modes. Further mod- 
elling is required to determine whether this is inconsis- 
tent with Kawaler's model. 



5 DISCUSSION 

While the model of line index oscillations we have pre- 
sented here is crude, it appears to explain the obser- 
vations quite well. We now must consider whether the 
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Figure 8. Balmer line line index with frequency 2742.85 (iRz, as 
a function of Balmer line, compared with a simple model (dashed 
line). The fit shown here gives <ST =~372K and Slogg =~0.047. 



Table 1. Changes in temperature and surface gravity for each 
frequency found in Paper II derived from Balmer line line indices 
for both 1999 and 2000 observations. 



Year 


Frequency 


AT eff 


Alogg 


fPaperll 




(/iHz) 


(K) 




(kms- 1 ) 


2000 


2742.85 


372±24 


0.047±0.006 


7.17 




2742.47 


196±13 


0.028±0.004 


4.45 




2102.83 


161±15 


0.018±0.004 


3.62 




2102.48 


392±24 


0.042±0.006 


8.47 




2101.57 


116± 7 


0.011±0.002 


3.40 




2075.72 


205± 9 


0.023±0.002 


4.27 




1985.75 


98±11 


0.013±0.003 


4.13 




1891.01 


67±11 


0.009±0.003 


1.99 


1999 


2742.63 


283±14 


0.051±0.004 


5.80 




2102.15 


560±17 


0.062±0.005 


11.20 




2075.29 


194± 4 


0.018±0.001 


2.41 




1890.98 


54±16 


0.003±0.004 


3.14 



changes in temperature and surface gravity we deter- 
mined arc reasonable. 



5.1 Temperature Changes 

We can use the Johnson B photometry presented in Pa- 
per I to test the plausibility of the temperature changes 
we have calculated. First, we need to convert these pho- 
tometric amplitudes to fractional bolometric luminosity 
changes. Kjeldsen & Bedding (1995) gave the following 
relationship between the fractional bolometric luminos- 
ity variation (5L/L)boi and the observed fractional lu- 
minosity variation (5L/L)\ 

T 



bol 



5L 

T 



A 

Abol 



(4) 



where 



A b 



ol 



he 



off 



(5) 
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Table 2. Change in radius determined from log g variations, com- 
pared with change in radius expected from velocities for a radial 
mode. 



Alogg 



Figure 9. AT c ff as a function of Alogg for the amplitudes shown 
in the top section of Table 1. The error ellipses were calculated 
using linear regression, and show Ict errors. The dashed line rep- 
resents a weighted least-squares fit to the data. 



Year 



Frequency 
(MHz) 



I AR\ ( AR\ 

V R /log 9 V R I 



(%) 



AR\ 
R > ve l 
(%) 



2000 


2742.85 


5.41 


0.44 




2742.47 


3.22 


0.27 




2102.83 


2.07 


0.29 




2102.48 


4.84 


0.68 




2101.57 


1.27 


0.27 




2075.72 


2.65 


0.35 




1985.75 


1.50 


0.35 




1891.01 


1.04 


0.18 


1999 


2742.63 


5.87 


0.36 




2102.15 


7.14 


0.90 




2075.29 


2.07 


0.20 




1890.98 


0.35 


0.28 



Equation 4 is a linearised expression, assuming the 
star is radiating as a blackbody. We have derived the 
exact form of this equation, where we also allow for the 
radius change in the star (in PG 1605+072, assuming 
radial modes, the radius change is 0.1-0.9%): 

(5L\ \{$L\ 8R~\ 

i(l-.-"-")t!?. (6) 
Abol n 

As an example, consider the 2102.15 y^Hz mode in 
our 1999 data (Paper I). The radius change is ~0.9% 
and the photometric amplitude is 3.2%, implying that 
(SL/L)b i ~10%. This leads to a temperature change 
of ^490K, which is roughly consistent with the value 
of 560 K shown at the bottom of Table 1. The tem- 
perature changes from photometry for the other modes 
from Table 1 are also about the same order of magni- 
tude as those determined from line index. We conclude 
from this that the temperature values we have derived 
are the correct order of magnitude. 

5.2 Surface Gravity Variations 

It is possible that we are measuring changes in the effec- 
tive surface gravity rather than the true surface gravity. 
The effective surface gravity is a combination of the true 
surface gravity and other forces acting on the stellar 
photosphere. In order to test whether we are measur- 
ing true surface gravity variations, we have calculated 
the radius change expected from the best-fit variation 
in surface gravity at each frequency. These changes are 
shown in Table 2, along with the radius changes cal- 
culated from velocity amplitudes, assuming all modes 
are radial. The difference between the radius change for 
a radial and non-radial mode is a scaling factor which 



is dependent on the inclination angle and limb darken- 
ing. For example, with no limb darkening and with the 
star pole on, an I = 1 , m = mode will have an ap- 
parent radius change ~30% larger than an I = mode 
with the same frequency, while an I = 2, m = mode 
will have a radius change ~10% smaller (Christensen- 
Dalsgaard 1994). Bedding ct al. (1996) presented an 
analysis for solar-like stars that also considered limb 
darkening. They found that the mode sensitivity of the 
line index of Balmer lines is similar to that of Balmcr- 
line velocity. 

All of the values of AR/R calculated from the log g 
changes are larger than those calculated from velocity 
by at least a factor of five. A more realistic interpreta- 
tion of the line index measurements requires sophisti- 
cated models. Phase dependent synthetic spectra need 
to be calculated that also account for limb darkening 
and rotation. As a prerequisite for such model spectra, 
the time dependent surface distribution of T e ff, log g 
and the velocity field have to be determined. Townsend 
(1997) has developed such a program (BRUCE) for the 
analysis of pulsations in rotating hot massive stars. Fal- 
ter (2001) recently wrote a code for modelling synthetic 
spectra of sdB stars from BRUCE output which shall be 
used to analyse the line index variations. A much larger 
data set has been acquired in May/ June 2002 by spec- 
troscopic and photometric multi-site campaign which 
became known as the Multi-site spectroscopic telescope 
(MSST, Heber et al. ). We will use the more detailed 
modelling when we shall analyse this much larger data 
set. 

However we may also have to check the validity 
of the hydrostatic approximation. Jeffery et al. (2001) 
encountered a similar problem in the analysis of the ex- 
treme helium star V652 Her; for this radially pulsating 
star the discrepancy between radii derived from surface 
gravity and spectrophotometry is about a factor of two. 
They found that at least for 10% of the pulsation cycle 
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of this star the assumption of hydrostatic equilibrium 
is not valid, due to shocks moving through the atmo- 
sphere. For PG 1605+072 the speed of sound is more 
than twice as large as for V652 Her due to its different 
chemical composition. A study of such effects are be- 
yond the scope of this paper and should be done after 
the more detailed hydrostatic spectral analysis outlined 
above has been carried out. 



6 CONCLUSIONS 

We have detected line index variations in the pulsating 
sdB star PG 1605+072. There is a strong dependence 
of line-index oscillation amplitude on Balmer line. We 
have developed a simple model, assuming hydrostatic 
equilibrium, where we use model spectra at as close to 
the same resolution and sampling of the observations as 
possible. Using this model, we have shown that despite 
not measuring true equivalent width variations, we can 
use the observables to infer effective temperature and 
effective surface gravity variations in PG 1605+072. 

Using more detailed models, we should be able to 
measure the AT and Alog g independently of photome- 
try and velocity measurements. With this information, 
we will have another set of amplitudes which we can use 
for mode identification of sdB stars. 

The authors would like to thank Mike Ireland for 
helpful discussions on error ellipses, and Steve Kawaler 
for continuing fruitful discussions on sdB evolution and 
interiors. This work was supported by an Australian 
Postgraduate Award (SJOT), the Australian Research 
Council, the Danish National Science Research Coun- 
cil through its Center for Ground-based Observational 
Astronomy, and the Danish National Research Founda- 
tion through its establishment of the Theoretical Astro- 
physics Center. 
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